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We calculate predictions from defect models of structure formation for both 
the matter and Cosmic Microwave Background (CMB) over all observable 
scales. Our results point to a serious problem reconciling the observed large- 
scale galaxy distribution with the COBE normalization, a result which is 
robust for a wide range of defect parameters. We conclude that standard 
scaling defect models are in conflict with the data, and show how attempts to 
resolve the problem by considering non-scaling defects would require radical 
departures from the standard scaling picture. 



PACS Numbers : 98.80.Cq, 95.35+d 

Defect models offer an elegant explana- 
tion of the origin of cosmic structure. The 
idea is that some distribution of defects— 
or more generally, field disorder — is pro- 
duced during a cosmic phase transition [|TJ. 
The defects then start a process of 'coars- 
ening' which continues through the present 
day, contributing a component to the mat- 
ter in the universe which evolves in a highly 
non-linear way. Cosmic strings, for ex- 
ample, move at relativistic velocities, peri- 
odically self-intersect and break off loops, 
which themselves eventually decay into grav- 
ity waves. Such processes can seed the onset 
of gravitational collapse in a universe which 
is initially perfectly homogeneous. 

In contrast with other models of cosmic 
structure formation, calculations for defects 
require the modelling of highly non-linear 
processes from very early times (e.g. the 
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time of Grand Unification), right up to the 
present day. During this period the universe 
increases by around 25 orders of magnitude 
in size. Only recently has it become prac- 
tical to solve the full Boltzmann equations 
for the matter and radiation perturbations 
in the presence of defect sources consistently 
modelled over such a length of time ||. Ac- 
curate large-scale numerical simulations are 
currently the best source of information on 
the details of the defect evolution, but even 
using state of the art technology it is still nec- 
essary to extrapolate with scaling arguments 
IJ to achieve anything like the required dy- 
namic range. Here we describe work that is 
not as closely linked to specific defect simula- 
tions and is therefore able to explore a wider 
range of possible defect models. Thus we can 
systematically investigate the robustness of 
the clash between defect models and obser- 



1 



vations. 

Our calculations use the fact that if only 
the power spectra are to be calculated, then 
one only needs two-point functions of the de- 
fect stress energy as input ||. Scaling argu- 
ments can then be used to increase the dy- 
namic range. This scaling behavior has been 
observed to some degree in numerical simula- 
tions and has become part of the stan- 
dard lore of defect evolution, although the 
extent to which it is valid over a factor of 
10 25 in cosmic expansion is not yet clear. 

For the current calculations this approach 
was incorporated into a version of CMBFAST 
H modified to include source stress energy 
for the scalar, tensor and vector contribu- 
tions, which are generic in defect based mod- 
els p iiT]|JTT[ ]. In order to do this, we model 
the components of the defect stress energy 
under a number of simple assumptions which 
maintain causality and conserve stress en- 
ergy. The source is approximated by a net- 
work of line-like segments with correlation 
length £77 at conformal time i] and velocity 
taken from a gaussian distribution with RMS 
v , truncated to prevent v > c. The number of 
lines is reduced causally, so as to maintain a 
constant density with respect to the horizon. 
This approach is similar to the model used in 
ref. ||12|| , which was shown to give two point 
functions in good agreement with string sim- 
ulations for certain stress energy components, 
but we have updated it to include all the 
components required and an improved decay 
mechanism flUf . 

Any active source which creates perturba- 
tions inside the horizon is likely to be inco- 
herent, leading to the absence or suppression 
of secondary Doppler peaks Hence, the 
form of unequal-time correlators (UETC) are 
also important. Our approach, which con- 
trasts with that used in ref. p|, is not to 
calculate the UETC directly, but to create 
an ensemble of source histories with the cor- 
rect two-point correlation statistics. Then in 



FIG. 1. The (COBE normalized) angular 
power spectrum of CMB anisotropics for the 
standard cosmic string model (solid) plotted 
with the current observational data, the stan- 
dard CDM curve (dotted). The two dashed 
curves give the partial contributions from two 
time windows to either side of z = 100 

order to calculate the ensemble average of 
the matter power spectrum and the CMB, 
one must use the Boltzmann code for each 
source history and average the resulting spec- 
tra. The results shown here used 200-400 re- 
alizations which gives very small statistical 
errors, but runs with just 40 realizations suf- 
ficed to establish the basic picture. 

In Fig. [I] we plot the angular power spec- 
trum of the CMB for what we shall call the 
standard string model (solid line). This uses 
string model parameters £ = 0.3 and v = 0.65 
as suggested by simulations, an assumption 
of perfect scaling from defect formation to the 
present day and a flat background cosmol- 
ogy with Q c = 0.95, Q h = 0.05 and h = 0.5 
where Hq = 100/ikmsec _1 Mpc _1 . Included 
also are the standard CDM model based on 
inflation (dot-dashed line) and all the current 
published data-points with error-bars based 
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on the assumption of gaussianity [15 



main features to note are the absence of any 
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discernible Doppler peak in the defect spec- 
trum and the apparent conflict with the data 
points. We have repeated these calculations 
for various different values of £ and v and 
also for sensible variations of the cosmological 
parameters h and fib- The spectrum is mod- 
ified by these variations, but none manage 
to increase the amount of power at angular 
scales with 1=200-400 by very much. Clearly 
the situation looks bad for defect models, al- 
though it is worth noting that the plotted er- 
ror bars are one-sigma, and deviations from 
the assumed gaussianity may require even 
larger error bars due to the small sky cov- 
erage. We expect the situation to be much 
clearer when the new CMB data arrives in 
the near future. 

Fig. |TJ also shows the partial results which 
come from integrating the defect contribu- 
tions over two time windows: Window 1 (z 
= 1300 to z = 100, during which rj increases 
by a factor of 5) gives the long-dashed curve 
and window 2 (z = 100 to z = 1.6, during 
which 7] increases by a factor of 7) gives the 
short-dashed curve. This information will be 
helpful in the subsequent discussion of non- 
scaling defect models. 

Using the same annotation as Fig. [I], 
Fig. |H shows the COBE normalized cold dark 
matter density perturbation power spectrum 
predicted for the standard string model and 
that for standard CDM, compared with the 
data [If]]. The contributions from the same 
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two time windows are also included, as in 
Fig. 1. Theory and data are often compared 
using <7g, the variance of the fractional mat- ^ioo ar e remote 
ter overdensity in a ball of radius 8/i -1 Mpc. 
For standard CDM a 8 = 1.2 for h = 0.5, 
while the value favoured by observations is 
erg = 0.5. If one were to compare with 
the string model at these scales, one calcu- 
lates 08 — 0.26 and hence the bias on these 
scales between the galaxy distribution, which 
is largely baryonic matter, and the Cold dark 
matter is b$ = cr 8 /cr° M ~ 2. This not unrea- 
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FIG. 2. The power spectrum of the dark 
matter perturbations for the same models and 
windows shown in Fig. 1 , plotted with the data. 

sonable value can even be slightly reduced by 
changing the string model and cosmological 
parameters. 

However, these comparisons ignore the 
fact that there is a woeful absence of power 
on larger scales. We quantify the conflict 
for low k by calculating the hypothetical bias 
b wo = ctioo / 0"k)o where <7i o is defined in anal- 
ogy to <t 8 , but for spheres of 100/t _1 Mpc and 
the favoured value (oioo = 3.5 x 10~ 2 ) is 
calculated for a smooth curve which gives a 
good fit to the data points. The standard 
string model has &ioo = 5.4 which cannot be 
improved substantially by any of the varia- 
tions already discussed. The chances of a real 
physical model having such a large value of 
nj], and there is no observa- 



tional evidence for a large &ioo j^8i|19| . We 
conclude that the standard string model is 
in conflict with the observations at an unac- 
ceptable level on scales around 100/i~ 1 Mpc. 
Variations relating to possible systematic un- 
certainties in simulations and our knowledge 
of the baryon density and Hubble constant 
cannot alleviate the discrepancy. We also 
note that these conclusions do not depend on 
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FIG. 3. The angular power spectrum of 
CMB anisotropics for the various non-scaling 
models discussed in the text. The three most 
extreme models (which have reasonable values 
of &ioo) have the highest peaks. Standard CDM 
is included for reference (dash-dot curve). 

the stringy nature of the model that we have 
used. For instance, the results are very sim- 
ilar if we impose a sharp sub-horizon cut-off 
on the source stress energy, mimicking be- 
haviour closer to that of cosmic textures. 

By far the most effective way of address- 
ing the large &ioo problem is to exploit the 
uncertainty in the overall scaling behavior of 
the string network. It is this behavior, after 
all, that relates the contributions from de- 
fects on different scales, and has simply been 
put into our calculation by hand. Although 
scaling has been observed to some degree in 
simulations, it is not completely clear that 
simple laws are valid over a wide dynamic 
range and events during the history of the 
universe may lead to deviations. For exam- 
ple, the radiation-matter transition is known 
at the very least to cause a shift; one could 
speculate that this is not yet well understood. 

We have extensively probed the possi- 
ble deviations from the standard picture, 
and found it very difficult to get around 
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FIG. 4. The power spectrum of the dark 
matter perturbations plotted with the data for 
the same models as Fig. || 

the large &i o problem. This can be under- 
stood by looking at the contributions from 
the two time windows illustrated in Figs. [I] 
and 0. The first window provides essen- 
tially all the contributions to the COBE nor- 
malization, while the second window pro- 
vides the dominant contribution to Oioo- The 
problem is that these two windows span a 
sufficiently narrow period in the defect his- 
tory that something dramatic must happen 
to the scaling behavior to shift their rela- 
tive contributions sufficiently. An extreme 
(and un-motivated) 'solution' which suggests 
itself is to simply turn off the string network 
at z = 100, hence preserving perturbations 
which contribute to Oioo, but removing the 
highest possible fraction of the contributions 
to COBE scales. The result of doing this is 
illustrated in Figs. [5] and f| (solid line), and 
manages to give b wo = 1.2. 

There are two simple types of deviations 
from scaling which may be more acceptable. 
The first, which has been observed to occur 
to some degree at the radiation-matter tran- 
sition P0]] , is just a step in the string den- 
sity from one value to another, occurring in a 



smooth way over some period of time. Such 
a deviation can be quantified by the ratio 
x = (r/ 1 / 2 ^o) rad /(^ 1/2 ^oo) mat (where x = 1 
gives "standard scaling"). For instance, if 
the radiation-matter transition gives rise to a 
difference in the amount of small scale struc- 
ture on the strings in the two eras, then \ 
could be interpreted as the ratio of the renor- 
malized string tensions. The second type of 
deviation we consider is a power law devia- 
tion from scaling quantified by a parameter 
a via #00 ~ r]~^ +a \ for which the density in 
strings p ~ l/r] 2+2a , with the choice a = 
corresponding to a standard scaling law. This 
may model the behaviour, for example, in an 
open universe or in one dominated by a cos- 



mological constant |21| . |22|j . 

In order to illustrate the problem, we 
show the results from four further models for 
the CMB in Fig. 3 and for the matter power 
spectrum in Fig. 4. The first two are mild de- 
viations from scaling which one might imag- 
ine are plausible: Model A (dotted curve, 
&ioo — 3.4) is a transition of x — 2, with 
the transition beginning at 8rj eq and ending 
at 10?7 eq where i] = r] eq is the time of equal 
matter-radiation. Model B (long-dash curve, 
&ioo = 2.9) is a power law deviation from scal- 
ing with a = 0.25. 

The other two examples are much more 
extreme; their virtue being that they can fit 
the data points in the matter power spec- 
trum at around 100/i _1 Mpc: Model C (long- 
short dashed curve, &i o = 1-0) is a transi- 
tion between the same times as for model A 
but with x — 10 and Model D (short-dashed 
curve, 6100 = 0.7) is a power law deviation 
from scaling with a = 0.75. Whilst models 
C and D fit the matter spectrum at scales of 
around 100/i~ 1 Mpc, they completely fail to 
fit smaller scale galaxy data, and give a large 
excess of small scale power in the CMB spec- 
trum. One might hope that changes to the 
ionization history and matter content of the 
universe could solve some of these problems. 



However, the fact remains that there is no 
evidence to suggest that such extreme devia- 
tions from scaling could occur in the standard 
defect scenario. 

We now relate our results to previous 
work on the subject. Several papers have 
discussed the bias in COBE normalized de- 
fect models. In ref. a serious bias problem 
was noted on scales up to 20/i -1 Mpc, but 
concerns remained that the simulations were 
not including all the relevant contributions 
(particularly to the density fluctuations) be- 
cause of their limited dynamic range. The 
compilation of refs. || and |23j in ref. pi"l , 
although looking very much like our Fig. 2, 
involved very different treatments of the de- 
fects at scales relevant to COBE vs feioo, and 
it is not clear that a straightforward compi- 
lation is valid. Our work avoids these uncer- 
tainties by solving the full Boltzmann equa- 
tions with a single source model to compute 
perturbations consistently on all scales. The 
scaling assumption translates into essentially 
'infinite' dynamic range. 

There has also been work more recently 
which is on equal footing in this respect ||. 
Conceptually, the main difference between 
this work and ours is that they extract the 
UETC's directly from numerical defect sim- 
ulations. Both groups scale the correlation 
functions to gain dynamic range. In addi- 
tion, ref. uses some sophisticated methods 
to work efficiently with the UETCs. Bear- 
ing in mind these differences, it should be 
noted that the two results look very similar. 
A strength of the simulation based approach 
is that the UETC's are associated with well 
defined defect scenarios. On the other hand, 
our approach is more flexible, allowing us to 
explore a wide range of UETC's in order to 
test the robustness of the results against vari- 
ation amongst defect models as well as pos- 
sible systematic uncertainties in the simula- 
tions. 

Even at a more technical level there is a 
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large degree of similarity between our results 
and those in ref. @. We find that on super- 
horizon scales the scalar, vector and tensor 
anisotropic stresses are in the simple ratio 



(i©i 2 ) = m 



2 ) : (|9jf) = 3 : 2 : 4 as 
imposed by causality and isotropy. Around 
/ = 10, we find that Cf : Cf : Cj approx- 
imately 3:1: 0.4. Our value of Cj : Cf 
is very close to that in ref. ||, while our 
Cf : CY is somewhat larger. The degree of 
similarity is striking given that we use a sim- 
ple model while large simulations were used 
in ref. @. We believe our larger Cf : Cf 
is due to the relatively large value of O o 
compared to Q s in our model. We have not 
directly compared our sources with those in 
|TT| (describing local strings), but our scalar 
component seems to be larger. These differ- 
ences (and those between |TT[ and [Q]) suggest 
that the relative strength of the scalar com- 
ponent may vary noticeably from one type 
of defect to another. We should note that 
even if there were only a scalar component 
the comparison with the current data would 
be very bleak; the vector and tensor compo- 
nents only make things worse. While it has 
been suggested that models with highly sup- 
pressed anisotropic stresses might achieve im- 
proved values for the COBE normalized bias 
|25| , no concrete defect model has been pro- 
posed which has this feature. 

Therefore, we conclude that standard 
scaling defect scenarios are ruled out on the 
basis of current data, and that this situation 
can only be remedied by extreme modifica- 
tions to the scaling law. 
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